Synthesis of superheavy elements beyond Z=118 
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We investigate the cold fission/fusion paths of superheavy nuclei within the two center shell model, 
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Z > 118. 
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The synthesis of superheavy elements beyond Z = 104, suggested by Flerov [l[, was predicted 
within the so-called fragmentation theory in Ref. by using the cold valleys in the potential energy 
surface between different combinations, giving the same compound nucleus. Soon it was shown in 
Refs. 0, 0| that the most favorable combinations with Z > 104 are connected with the so-called Pb 
potential valley, i.e. the same valley of the heavy cluster emission 

Due to the double magicity of ''^Ca, similar with ^°*Pb, in Ref. [3l it was proposed ''^Ca as a 
projectile on various transuranium targets. Indeed, the production of many superheavy elements with 
Z < 118 (corresponding to the last stable element Cf) during last three decades was mainly based on 
this idea [ill . 

The formation of superheavy compound systems by fusion was intensivelly explored (lol . [ill , [l^ . 
On the other side the investigation of experimental data concerning fusion and fission of superheavy 
nuclei with Z = 112, 114, 116, together with data on survival probability of these nuclei in evaporation 
channels with 3-4 neutrons, revealed the fact that the fission barriers are quite high, leading to a 
relative high stability of such systems [l3| . 

The main tool to investigate such nuclei is almost exclusively based upon the investigation of a- 
decay chains. In the last decade several papers were devoted to the calculation of a-decay half-lives in 
this region [l3|. AH these approaches can be considered as phenomenological ones, based essentially 
on the Gamow a-nucleus potential picture p^ . The recent microscopic estimate of the a-particle 
preformation factor, by using shell model single particle orbitals, performed in Ref. [lH, showed 
that the strong change of the Q-value along neutron chains can by explained only by supposing the 
existence of an a-cluster component in heavy and superheavy emitters. 

The a-particle emission is connected with the " lightest" side of the cold valley on the fragmentation 
potential surface. On the other hand the "heaviest" side of the cold valley is given by the cold fission 
process, i.e. the emission of two fragments with similar masses in their ground states. Between theses 
limits there is a broad region of cold heavy cluster decays. The aim of this work is to evidence cold 
fission valleys, which can be good candidates for the production of superheavy elements with Z > 118, 
by using the inverse, fusion process. We extend the analysis performed in Ref. [TtI within a simple 
model, to a more reliable microscopic approach to estimate the fission/fusion barrier, given by a new 
version of the Super Asymmetric Two Center Shell Model [l^l ■ This version solves a Woods-Saxon 
potential in terms of the two center prescriptions and provides two additional degrees of freedom, 
that is, the deformations of the fragments. The deformation energy of the di-nuclear system is the sum 
between the liquid drop energy and the shells effects, including pairing corrections. The macroscopic 
energy is obtained within the framework of the Yukawa-plus-exponcntial model extended to binary 
systems with different charge densities [20|. Strutinsky correction prescriptions 21 1 were computed on 
the basis of a new version of the two center model with Woods-Saxon deformed potentials of fragments. 
We considered only cold fission/fusion process. Consequently the deformations of the initial and final 
nuclei are given by their ground state values of Ref. [22] . 

The penetrability, coresponding to some binary partition, defines the isotopic yield and it is char- 
acterized by the difference between the nuclear plus Coulomb potential and the Q-value. For a given 
initial nucleus (Z, A) this quantity, called drivind potential, depends upon the charge, mass numbers 
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of a given fragment (we will consider the second one) and the inter-fragment distance. For a fixed 
combination A — Ai + A2 the driving potential has a minimum at the charge equilibration point Z2, 
which we will not mention in the following, i.e. 

V{A2,R) = Vn{A2,R) + Vc{A2,R + B{Zi,Ai) + B{Z2,A2) , (1) 

where Vn{A2,R) is the nuclear and Vc(^2,^) Coulomb inter-fragment potential. Here the binding 
energy of the initial nucleus is not considered, because it has the same value for all binary partitions. 
For deformed nuclei, due to the fact that the largest emission probability corresponds to the lowest 
barrier, this potential decreases in the direction of the largest fragment radius. We mention that a very 
convincing theoretical evidence that the cold fission has a sub-barrier character and depends upon 
deformation parameters was the calculation of penetration factors, using a double the double folding 
inter-fragment potential, with M3Y plus Coulomb nucleon-nucleon forces. This simple estimate was 
able to reproduce the gross features of the binary cold fragmentation isotopic yields from ^^^Cf [23[ . 
Here it was shown the important role not only of quadrupole, but also of hexadecapole deformation 
parameters on emission probability. 

Our first step is to estimate the driving potential ([1]) for nuclei beyond Z=118. First of all we have 
checked the validity of the two center shell model, by reproducing qualitatively cold isotopic yields 
from ^^^Cf [2^. Then, we computed the driving potential for the hypothetical compound nucleus 
■^20^ ■ This nucleus has a larger amount of neutrons than the measured superheavy combinations with 
Z=116 and 118 24]. In Fig. [T]we plotted the two dimensional potential surface V{A2,R) versus the 
fragment mass number of the second partner and the inter-fragment distance R. It is clear that the 
cold fission/fusion process is very much hindered in the region of the large maximum between A2 = A 
and = 40. 




FIG. 1: The two dimentional potential V{A2,R) versus the fragment mass number and the inter-fragment 
distance R, for the hypothetical compound nucleus 120^. 



This is especialy clear from Fig. [SJ where we plotted the maximum value of the potential surface 
in Fig. [T]with respect to the inter- fragment radius R, as a function of the fragment mass number. 

The first relevant minimum of the potential surface in the region A2 > 40 corresponds to the already 
mentioned double magic nucleus "^^Ca, but unfortunately its binary partner ^^^Fm is unstable, as can 
be seen from Table 1. 

In order to search for reliable binary candidates, producing the above mentioned hypothetical com- 
pound nucleus, it is necessary to investigate the penetration factor. This quantity can be estimated. 
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FIG. 2: The maximum value of the potential surface in Fig. Q] with respect to the inter-fragment radius R, 
versus the fragment mass number. 



as usually, by using the semiclassical integral 

Pa, = exp I -2 [ViA2, R) - B^Z, A)]dR^ , (2) 

between internal and external turning points. In Fig. [3] is given the penetrability ratio R = PA2/Pa, 
corresponding to the potential surface in Fig. [T] Several maxima, comparable or larger than the 
mentioned combination ^*Ca+^^^Fm, are present in this figure. They correspond to Cr, Fe, Ni, Zn, 
Ge, Se projectiles. The most promising are the projectiles beyond Zn isotopes, where the penetrability 
increases by nine orders of magnitude from ''^Zn to '''^Zn, as can also be seen from Table 1. The region 
of maximal values corresponds to Sr+Pb combinations, i.e. to the already mentioned Pb valley. 

On the other hand the target-projectile combinations should be relative stable. In Table 1 are 
given the fragment half-lives of these binary combinations, with A2 > 40, corresponding to local 
maxima of the penetrability. Unfortunately the half-lives of ^^'^^Zn isoptopes and their partners are 
fast decreasing. 
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FIG. 3: The penetrabilities corresponding to the potential surface m Fig. Q] versus the mass number of the first 
fragment. 



Table 1 

Mass numbers and half-lives of binary partners for the hypothetical compound nucleus 120^- In the 
last column is given the penetrability ratio R = P^^ /Pa . 



Ai Y 
Zi ^1 


Ti 


A2 Y 
Z,X2 


T2 


R 


48 Pa 
2o>^d 


1019 y 


252Fm 


25.4 h 


2.1 


10- 


26 




13.9 s 




30 m 


3.1 


10" 


25 


52 rp- 

22 -LI 


1.7 s 


248 pf 
98 


333.5 d 


9.5 


10- 


26 




5.94 m 


244 Pm 
96 


18.1 y 


3.7 


10' 


24 


58 


7 s 


242 
96 


162.8 d 


1.7 


10- 


23 


60 


490 ms 


240 
96 


27 d 


5.5 


10- 


26 




160 ms 


238p™ 
96 


2.4 h 


4.8 


10- 


21 




2 s 


236p„ 

94 


2.9 y 


6.3 


10" 


30 


iNi 


54 h 


234TT 
92 ^ 


2.5 10^ y 


2.5 


10- 


25 


iNi 


29 s 


232TT 
92 ^ 


68.9 y 


5.4 


10- 


23 


2'iNi 


1 /iS 


230TT 
92 ^ 


20.8 d 


8.1 


10- 


22 


30^1^ 


46.5 h 


228t}i 
90 ^'^ 


1.9 y 


3.7 


10- 


20 


30 


95.6 s 


226rpi 

90 -Lil 


30.6 m 


5.7 


10- 


11 



Ax Y 




A2 Y 


T2 


R 




1025 y 


88 


3.7 d 


8.1 


10- 


-11 


iGe 


88 m 


222n 

88 


38.0 s 


7.1 


10" 


-08 


i2Se 


stable 


220t3„ 

86 


55.6 s 


1.4 


10- 


-11 


82 
34^6 


1019 y 


218r) 

86 


35 ms 


5.2 


10- 


-12 


84 n 

34oe 


3.1 m 


216r) 

86 


45 ^s 


4.1 


10- 


-15 


86 Op 
34^6 


15.3 s 


214r)„ 

86 


270 ns 


1.0 


10- 


-12 


88Kr 


2.84 h 


212p 
84 ^0 


299 ns 


3.9 


10- 


-14 


36^1^ 


32.3 s 


210p 
84 ^0 


138.4 d 


7.4 


10- 


-08 


92 Ty- 
36^1^ 


1.84 s 


ifPo 


2.9 y 


8.1 


10- 


-08 


94 0-. 
38 '-'1 


75.3 s 




stable 


1.4 


10- 


-07 


96 
38 '-'1 


1 s 


i4pb 


stable 


2.2 


10- 


-04 


98 0^ 
38 


653 ms 


i^Pb 


5.2 10^ y 


1.1 


10- 


-02 



In order to compare the data in Table 1 with those corresponding to some experimentally detected 
superheavy elements, in Table 2 we give penetrability ratios of three measured fusion reactions. 

We also compared the penetrabilities in Table 1 with the corresponding quantities of some close 
elements which have been experimentally reported by using 48Ca as projectile. The reactions given in 
the first and second lines of the Table 2 have penetrability ratios around R « 10~24 The corresponding 
values of the first two lines in Table 1 are not far from this number. 



Table 2 

Mass numbers and half-lives of binary partners for the experimentally detected compound nuclei '^X . 
In the last columns is given the penetrability ratio R = Pa2 I Pa cind the quoted reference. 



Ai Y 

z, 


Ti 


A2 Y 
Z2X2 


T2 


iX 


R 


Ref. 


48 Ta 

20 '-^a 

48 Ta 

20 '-^a 

70 7n 
30 ■''11 


1019 y 

1019 y 

stable 


244 p 
96 

246 rif 

98 

ifPb 


18.1 y 
35.7 h 
stable 


292 Y 
116^ 
294 Y 
118^ 

278 Y 
112^ 


1.7 10-24 
1.1 10-24 

5.0 10-11 


[25] 
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The penetrability ratio, corresponding to the third reaction ^°*Pb+™Zn^f]'|X [2^, is R = 
Pa^/Pqi ~ 5 10~^^. PracticaUy all combination beyond ''''Zn, giving the hypothetical element f^o-^i 
have the penetrability ratios larger than this value. In Table 1 we remark that in the combination 
''^Ge+^^^Ra both partners have enough large half-lives to be used in fusion experiments. 

Concluding, wc computed the potential energy surface for different binary combinations, giving the 
superheavy compound nucleus 120X, by using the Two Center Shell Model. We evidenced binary part- 
ners whose penetration factors are comparable with similar combinations, producing close superheavy 
elements which were experimentally measured. An extensive analysis on fission/fusion channels for all 
possible isotopes in this region of superheavy elements, by using this performant version of the Two 
Center Shell Model, is under way. 
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